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We study the strong localization of atomic matter waves in a disordered potential created by atoms 
pinned at the nodes of a lattice, for both three-dimensional (3D) and two-dimensional (2D) systems. 
The localization length of the matter wave, the density of localized states, and the occurrence of 
energy mobility edges (for the 3D system), are numerically investigated as a function of the effective 
scattering length between the atomic matter wave and the pinned atoms. Both positive and negative 
matter wave energies are explored. Interesting features of the density of states are discovered at 
negative energies, where maxima in the density of bound states for the system can be interpreted in 
terms of bound states of a matter wave atom with a few pinned atomic scatterers. In 3D we found 
evidence of up to three mobility edges, one at positive energies, and two at negative energies, the 
latter corresponding to transitions between extended and localized bound states. In 2D, no mobility 
edge is found, and a rapid exponential-like increase of the localization length is observed at high 
energy. 
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I. INTRODUCTION 



The propagation of waves in disordered systems is a 
rich physical phenomenon, the object of enduring re- 
search interest. Its complexity is due to the fact that 
the scattering of a wave by a random potential depends 
on several features: energy and type of the wave, internal 
and external degrees of freedom of the scattering poten- 
tial, dimensionality and symmetry of the physical system, 
and possible presence of interaction among the propagat- 
ing waves. Such variety is the reason for a wide experi- 
mental and theoretical studies of the diffusion of several 
kind of classical and quantum waves by many kinds of dis- 
ordered potentials. This field of research was started by 
P.W. Anderson [l[, who predicted in 3D the localization 
of a quantum particle experiencing short range hoping 
between discrete sites, when the on-site energies are suf- 
ficiently random. Several theories have been developed, 
leading to a general consensus on the fact that a wave is 
localized in ID and 2D infinite disordered systems, inde- 
pendently of its energy and of the strength of disorder. 
In 3D systems a metal-insulator phase transition induced 
by sufficiently strong disorder may take place: there ex- 
ists one (or more) critical energy E c (called the mobility 
edge) which separates two energy regions. Waves with an 
energy on one side of E c are spatially localized (transport 
is absent), while waves with an energy on the other side 
are extended over the entire space (transport is diffusive) . 
The localized waves are characterized by an amplitude 
which decreases exponentially in space at large distances 
from a central region, defining a typical length, called the 



localization length £. The richness of the phase transition 
appears in the critical region around the critical energy 
E c : here the physics is supposed to be universal, de- 
pending only on the symmetries of the system, and not 
explicitly on the kind of wave or disordered potential. 
This is due to the fact that, in the critical region, the 
localization length diverges with a power law behavior 
£ oc \E — E c \~ v , where the critical exponent v character- 
izes the universality class of the phase transition ■ 

The localization of waves appears as a very rich phe- 
nomenon, but it is also very delicate and complex, often 
difficult to observe experimentally due to parasitic effects 
(absorption of the wave, finite size effects, ...). A large 
scientific literature exists on wave localization for sev- 
eral physical systems. Here we deal with the localization 
of ultra-cold atomic matter waves. Indeed, ultra-cold 
atomic gases offer a unique system in terms of isolation 
from environment, cleanness of the system, realization of 
the system in several spatial dimensions (ID, 2D, 3D), 
development of several direct detection techniques, and 
control of the interaction strength. One can even imagine 
localizing different kinds of waves (matter or light) in dif- 
ferent kinds of disordered potentials (created by light or 
matter). Recently, several experimental investigations on 
the localization of atomic matter waves have been per- 
formed: in a ID genuine disordered potential made by 
a laser speckle in a ID bi-chromatic optical lattice 
0, and in an atomic quasi-periodic kicked rotor where 
a localization in momentum space has been reported Q 
(true disorder is absent, but the system, which is ID in 
real space, can be mapped into an Anderson model with 
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several effective dimensions in momentum space, giving 
the possibility to extract the critical exponents of a 3D 
Anderson transition). 

A different way to realize a disordered potential for 
atomic matter waves (atoms of species A), using atoms 
(another species B) pinned at random positions at the 
nodes of an optical lattice, was proposed in This 
proposal has potentially several advantages over the laser 
speckle route. On the experimental side, it may be real- 
ized in principle as easily in ID, in 2D, and in 3D, and 
also since an optimization of the matter wave interac- 
tion with each individual scatterer may be performed, 
one may hope to reach very short mean free paths and 
localization lengths. On the theoretical side, it allows an 
exact numerical study for a large number of scatterers 
- as many as in typical experiments as we shall see. A 
quantitative study of this model was done for a ID sys- 
tem 0], predicting minimal localization lengths < 10/im. 
A first investigation of the model for 3D systems was per- 
formed in @, revealing the existence of a large density 
of localized states with a very short localization length 
(a few microns) for positive energies. Recently, similar 
models have been considered in the case where the mat- 
ter wave experiences both a per iodic potential and inter- 
actions among the A atoms |10| . 

In this paper, we study more extensively the behavior 
of the matter wave localization for this model, for both 
3D and 2D systems, for both positive and negative ener- 
gies, and for a broad range of values of the matter wave- 
pinned atom interaction strength. We explicitly study 
the features of the localization length and of the density 
of localized states, the occurrence of mobility edges, and 
the effect of the presence of an underlying periodic lattice. 
We find that, both in 2D and 3D, a strongly resonant in- 
teraction between the matter wave and each pinned atom 
is required to obtain localized states with small localiza- 
tion lengths. On the experimental side, this requires a 
tuning of this interaction by use of a Feshbach resonance. 
First steps in this direction were recently taken: confine- 
ment induced resonances between almost free A = 87 Rb 
atoms and B = 41 K atoms ti ghtl y trapped in a ID laser 
standing wave were observed [ll| . 

The article is organized as follows. The first part deals 
with the 3D case: We present the model in section [Til we 
calculate the localization length in section|nI]by studying 
the response of the system to a source of atomic matter 
waves A in the medium consisting of pinned B atoms, 
we calculate and discuss the density of localized states in 
section [TVT, in the form of resonances at positive energies 
and bound states at negative energies. The second part 
of the article contains the same analysis for the 2D case, 
see section [V] We conclude in section l*VTl 



II. PHYSICAL SYSTEM AND MODEL 

We consider two atomic species, A and B. The B 
atoms are tightly trapped at the nodes of a cubic opti- 




FIG. 1: (Color online) For a source placed at the origin ro = 
0, coarse grained histogram giving (ln|D|) as a function of 
the distance from the source, averaged over 100 realizations 
of disorder. The lattice for the B atoms is cubic with a lattice 
constant d, and a filling factor p occ — 1/10. The energy of 
the emitting source is E = h 2 k 2 /2m with kd = 0.94. The 
effective scattering length is a c g = 0.7d. (a) Coarse grained 
histogram for 5 different radii R of the sphere containing the 
scatterers: R/d = 40 (black), R/d = 50 (red), R/d = 70 
(green), R/d = 90 (blue), R/d = 110 (violet), (b) For R/d = 
110, coarse grained histogram (black circles) and its fit by the 
functional form Eq. ()10[) (black lower solid line), resulting in 
nd — 0.086 and a = 2.5; the same histogram (red circles) and 
the same fit (red upper solid line) with Di multiplied by rf . 
The fit was performed over the interval 10 < r/d < 90. Note 
that the largest considered value of R/d = 110 corresponds 
to a number of scatterers N ~ 5.6 x 10 . 



cal lattice of lattice spacing d, in a regime where their 
tunneling among the neighboring lattice sites is negligi- 
ble over the duration of the experiment. These B atoms 
are prepared in the vibrational ground state of the local 
microtrap, and arc randomly distributed among the lat- 
tice sites, with a uniform occupation probability p occ in- 
side a sphere of radius R; their locations are independent 
random variables, except for the constraint that there is 
nowhere more than one B atom per lattice site [i~2| . 

The A atoms form the matter wave to be strongly lo- 
calized. They are assumed to move freely in space (in par- 
ticular they are insensitive to the optical lattice), except 
that they scatter on the B atoms. This scattering is as- 
sumed to be elastic, under the condition that the kinetic 
energy of the A atoms is much smaller than the quan- 
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turn of oscillation of the trapped B atoms, <C Huj osc , 
where m and fi-fc are the mass and momentum of an A 
atom. This scattering is also assumed to be in the zero- 
range regime fc|r ff| <C 1 where r e ff is the effective range 
for the scattering of an A atom on a trapped B atom. 
Since r c s — a,ho when the modulus of the 3D free space 
A — B scattering length a is much larger than the har- 
monic oscillator length aho of the B atoms (see 0), we 
assume that ka^ <C 1. In this zero-range regime, the B 
atoms may be considered as point-like scatterers, and the 
A — B scattering is characterized by the effective scat- 
tering length a e s fl3j | . In practice, at positive energy 
E = h 2 k 2 /2m > of the matter wave atoms A, we shall 
consider at most an energy of the order of Eq = h 2 /md 2 , 
so that k < 1/d, and the zero-range condition becomes 
d 3> aho, a condition well satisfied in a deep optical lat- 
tice. At negative energies E < 0, an AB dimer of en- 
ergy E = -Edim = — / (2 ma es) exists for a e g > 0, such 
that the zero-range condition becomes a e s ^> aho- Re- 
markably, using a Feshbach resonance technique to adjust 



the free-space A — B scattering length a, one can real- 
ize confinement-induced resonances leading to arbitrarily 
large values of a c s @, so that the condition a c ff 3> &ho 
may be realized. The B atoms may thus effectively con- 
stitute a static and strong disordered potential for the 
matter wave. 

The problem is thus modeled as follows: The mat- 
ter wave Hamiltonian is that of the free A atom, % = 
— ^-A r , where A r is the Laplace operator, and the A—B 
interaction is replaced by the following contact conditions 
for the A atom wave- function ip(r): there exist complex 
numbers Di such that 

tf(r) = A D r - r *r 1 " acff 1 ] + 0(\r - r,|) (1) 

in the vicinity of each B scatterer location iv The fac- 
tor m/(2irh 2 ) is introduced for convenience. The Bcthc- 
Pcicrls contact condition (TTJ) is equivalent to the pseudo- 
potential as used in Q. 
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FIG. 2: (Color online) For the 3D system: Comparison of the two methods of calculating the Lyapunov exponent k = 
As in Fig[TJ the filling factor is p occ = 1/10, and the effective scattering length is a a g = Q.7d. (a) Lyapunov exponent k for a 
fixed energy of the emitting source, E — h 2 k 2 /2m with kd — 0.94. Red squares: Lyapunov exponent from the first method, see 
Eq. (|10p . for several values of d/R; each point is the result of a fitting on the interval lOd < r < R — 20d. The second method, 
see Eqs. (|13|l . (j 14 j> . and (|15p . leads to a Lyapunov exponent which is the extrapolation at d/R — of the black circles; the green 
straight solid line is the linear extrapolation, giving nd = 0.089 ± 0.005 at d/R = 0; the three horizontal dotted lines represent 
the value of k and its confidence interval. The two methods are essentially compatible if one considers the scatter in the data 
for the first method and the error bars for the second one. (b) For the two methods, Lyapunov exponent k as a function of the 
energy of the emitting source E — h 2 k 2 /2m. The lattice spacing d of the cubic lattice, and the energy Eq = h 2 /(md 2 ), are used 
as units. Note that Eq = 2E Tec /ir 2 , where _E rec is the atomic recoil energy after absorption of a lattice photon. Here the filling 
factor is p OC c — 1/10, and the effective scattering length is a e s = 0.7d. The solid lines and the violet star are obtained with the 
first method [see Eq. {TO}]. The three solid lines correspond to sphere radius R — 50d (black line, triangles), R — 60d (red line, 
squares), R — 70d (green line, circles), and an average over 500 realizations of disorder. The violet star corresponds to a bigger 
sphere with radius R = llOd, and 100 realizations of disorder. The blue diamond symbols with error bars, are obtained with 
the second method [see Eqs. (|13l) . (|14[) . and (|15[) ] where the extrapolation has been performed using d/R — 1/50, d/R = 1/60, 
d/R = 1/70, and 500 realizations of disorder. The second method provides results essentially in agreement with those of the 
first method. Note the quite small values of the Lyapunov exponents (large localization lengths). 

I 



Here, the disordered potential has a finite extension, so that all the positive energy eigenstates are extended 
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states belonging to a continuum. However, the localized 
states that would exist for a strictly infinite extension 
disorder have precursors in the form of sharp resonances 
with a width tending exponentially to zero with the dis- 
order extension dig. At negative energies the matter 
wave is bound inside the gas of scatterers, and the cor- 
responding bound states can be either extended or local- 
ized. The appropriate tool to find these resonances and 
bound states is the matter wave Green's function for an 
energy E, for a given realization of the disorder, 



G(r,r') 



1 



E + i0+ - H 



(2) 



and its analytical continuation to complex energies in the 
lower half-plane. Here we have used the standard nota- 
tion G(E + iO+) = lim,_> + G(E+ir)). Note that G(r,r') 
is subject to the same contact conditions as ip in ((T|), 
hence the subscript c.c. in Eq.@. This Green's function 
G actually plays a major role in the theory of transport 
phenomena @: The large |r — r'| behaviors of (G(r,r')) 
and (In |G(r, r')|) provide the scattering mean free path 
and the localization length, respectively, where (. . .) rep- 
resents the average over all realizations of disorder, that 
is over the B locations. 

A remarkable feature of the point-like scatterers is that 
the Green's function may be obtained by the solution of a 
N x N complex linear system, where N is the number of 
B scatterers. This allows an exact numerical calculation, 
up to N « 5 X 10 5 in this work, a value of the order of 
typical experiments with atoms in optical lattices. 

This remarkable feature may be derived as follows. 
One starts with the fact that the Green's function 
G(r, ro) obeys a Schrodinger equation with a point-like 
source term 5(r — ro) of matter waves at position ro. 
The diverging terms of the contact conditions ([1} for the 
Green's function give rise to secondary point-like sources 
of amplitudes Di at the scatterers' positions r^, by virtue 
of the usual relation A r |r — r^ -1 = — 47r<5(r — rj). The 
resulting wave equation is thus 



E + i0 + + — A r ) G(r, r ) = 5{r - r ) 



JV 



J2DiS(r-ri). (3) 



This may be integrated using the Green's function <?o( r ) 
for H in the absence of scatterers. We set 
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E = 



h 2 k 

2m 



(4) 



For E > 0, we impose k > 0, and 



5o (r) 



2irh 2 



(5) 



For E < 0, we take k = iq, with q > in Eq. ©. The 
resulting solution of Eq. (|3|) in thus 



N 



G(r,r ) =£fo(r-r ) + ^Affo(r-ri)- 



(6) 



The secondary sources amplitudes Di are then deter- 
mined by imposing on ([B]) the contact conditions (JTJ) at 
the order O(l). That is, for the non-diverging term l/a e ff 



N 



2TTH 2 



fldfa-ro), Vie {!,...,#} (7) 



where we have introduced the N x N matrix 



My = 



2nh 2 . . 
9o(*i -rj) if i ? 3, 

(8) 

ik + a~g if i = j. 



Eq. ([7]) constitutes the aforementioned N x N linear sys- 
tem, and its formal solution gives Q 



G(r,r ) = 5o(r-r ) 



2nh' 



N 



-ro 



(9) 

The matrix M will play a crucial role in what follows. 
Indeed, we will relate the localized eigenmodes of the 
matter wave to the eigenvectors of the matrix M. This 
matrix shows a purely off-diagonal disordered coupling 
between different scatterers, which is long-range at posi- 
tive energies, and short-range at negative energies. This 
kind of disorder does not coincide with the one originally 
introduced by Anderson, characterized by diagonal dis- 
order and short range couplings Finally, it is worth 
noticing that this model can be exactly mapped to the 
case of a scalar light wave scattered by a disordered en- 
semble of two-level atoms, where the same matrix M 
appears fl5| . 



III. LOCALIZATION LENGTH source is expected to be exponentially suppressed (lo- 

calized) at large distances, if its energy is smaller than 
A. Defining and calculating the localization length a critical energy E c , the so-called mobility edge. This 

In an infinitely extended disordered three-dimensional 
system, a matter wave of positive energy emitted by a 
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FIG. 3: For the 3D system: Coarse grained histogram (ln|7J]} 
as a function of the distance from the source, for two differ- 
ent kinds of disorder. The first kind of disorder (solid line) 
is obtained by a random filling of the cubic lattice inside the 
sphere of radius R = 70d. A second kind of disorder (dashed 
line) is obtained first by a random filling of the cubic lattice 
inside only one sector (1/8) of the sphere, x > 0, y > 0, z > 0, 
and then by filling the remaining seven sectors by reflecting 
the positions in the first sector with respect to planes x — 0, 
y = 0, z = 0. The average is taken over 50 realizations of 
disorder. The filling factor is p occ — 1/10, the effective scat- 
tering length is a e a = d and the energy of the emitting source 
is E = h 2 k 2 /2m with kd = 0.1. The circles represent, for the 
two kinds of disorder, the fit of the histogram on the interval 
lOd < r < 55d to extract the Lyapunov exponent k with the 
first method [see Eq. iflOjl]. We obtain [nd = 0.31, a = 2.2), 
and (Kd = 0.31, a — 3.6), for the first and second kind of 
disorder, respectively. The Lyapunov exponents appear to be 
identical, so that one can use the second kind of disorder to 
calculate k with a substantial gain in the computational effort 
(a factor 64 on the memory size and a factor 512 in the CPU 
time) . 



corresponds to absence of matter wave transport. The 
localization length £ is the length scale associated with 
this exponential decay. It is an average quantity, which 
has to be calculated by taking the mean over all possible 
realizations of the disorder. Waves emitted by a source of 
energy E larger than E c are instead expected to propa- 
gate diffusively in the system, with an intensity decaying 
as the inverse of the distance. 

As we shall see, mobility edges can also be present for 
negative energies. Indeed, a single matter wave A atom 
exhibits bound states with a few B scatterers (dimers 
AB, trimers AB-2, tetramers AB3, etc.), that have a 
non-zero hopping amplitude among different scatterers. 
Then, the evanescent matter wave emitted by a source at 
negative energy may populate such bound states. This 
may lead to waves of dimers, trimers, etc., bounded in 
the volume occupied by the scatterers, but which can be 
either extended or localized within that volume. 

To obtain these properties we thus calculate the 
Green's function G(r, tq) which gives the r-dependent 



matter wave amplitude resulting from the source in Tq. 
In practice we solve numerically the linear system Eq. ([7]) 
and we extract the localization length £ for the ampli- 
tudes Di in two different ways. 

First method: The spatial localization of the Green's 
function is reflected in a localization of the secondary 
source amplitudes Di. We thus calculate the average 
over disorder (ln|Dj|), and construct a coarse grained 
histogram of the data set (|rj — ro|, (In |£)j|}), that we fit 
with the functional form 



r-r 



In 



C- 



r- r o 



(10) 



where C, a and the Lyapunov exponent k are the three 
free parameters. The localization length is then £ = 1/k. 
The fit is performed over the secondary source amplitudes 
Di with positions such that r m - m < |r» — ro| < 7" max , in 
order to exclude the near-field contribution of the source, 
and to minimize the effects of the boundaries of the dis- 
order. Both effects appear in FigQJ,, which reveals a 
boundary layer of « 15 lattice spacings. Apart from these 
boundary effects, the results for increasing radii R of the 
sphere containing the scatterers are in good agreement, 
and reveal that the decay of Di is not simply exponential. 
On the contrary, the inclusion of the power-law factor in 
Eq. (flU)) provides an excellent fit to the numerical data, 
as shown in FigfT}). For the system sizes that we are 
able to treat numerically, the failure to include a as a 
fitting parameter would lead to unreliable values of k de- 
pendent upon the system size. Finally, we note that the 
functional form Eq. (flU)) includes both the localized and 
the diffusive sides of the phase diagram, with £ = +00 in 
the diffusive regions. 

Second method: Inspired by the usual definition of the 
Lyapunov exponent in one-dimensional disordered sys- 
tems, we define a direction dependent transmission coef- 
ficient t(n) for the field emitted by the source, for a given 
realization of disorder: 



G(r,r ) ~ t(n)g (T-r ), 

r— »+oo 



(ii) 



with the unit vector n = r/r. Here / ~ g means f/g—t 1. 
This expression, compared with Eq. (JB]) for r — > +00, 
provides the exact relation 



N 



t(n) = l + ^Ae" fc(r '" ro) 



(12) 



i=i 



In the one-dimensional case, the quantity to consider is 
the logarithm of the transmission coefficient, rather than 
the coefficient; the former is indeed a self-averaging quan- 
tity contrary to the latter [l6| . We thus define the direc- 
tion dependent three-dimensional Lyapunov exponent as 



K (n) 



_ Hm H*(n)l) for £ >0i (13) 
fl-^+oo R 

q lim Qa|t(n)l) i0lE<Qi (14) 

i?,->+oo R 
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FIG. 4: (Color online) For the 3D system: Localization length £ [(a)], and Lyapunov exponents « = l/£ [(b)], as a function 
of the energy E of the emitting source for different values of the effective scattering length a e fl?: a ff = 0.1a! (black solid line), 
ttcff = 0.2a! (red solid line) , a e s = 0.7a! (green solid line), a e ff = d (blue solid line), a e g — 1.3a! (violet solid line). The mean 
field calculation [see Eq. (|f 8[) ] is also shown for a e a = 0.1a! (black dashed line), and for a e s = 0.2a! (red dashed line). In (a), 
for E/Eo = —1.2, one has a e g/d = 1.3, 1, 0.7 from bottom to top, and for E/Eo = 0.4, one has a e g/d = 1, 0.7, 1.3 from bottom 
to top. The lattice spacing d of the cubic lattice, and the energy Eo = ft 2 /(md 2 ), are used as units. Results are obtained using 
the first method [see Eq. (|10|) ]. with 500 realizations of the disorder for E > and 100 realizations for E < 0. The filling factor 
is Pocc = 1/1, and the sphere radius is R = 70a!, which leads to a mean number of scatterers (N) ~ 1.4 x 10 J . Comparing the 
results for E > with 100 and 500 realizations, we estimate that the error on « is « 10%, and even smaller for the lowest 
values of £. 



where R is the radius of the sphere containing the point- 
like scatterers and the average is taken over the disorder 
[It}. The term q in fTJJ has been introduced since, for 
E < 0, the field emitted by the source go(r — ro) taken 
as a reference field in Eq. fTTl) decays as e~ qr for r —> oo. 

We have found numerically that At(n) only weakly de- 
pends on the direction so that we may define the local- 
ization length in terms of its average over solid angle: 

«=^/^(n). (15) 

In practice, we calculate the limit in Eq. (flU]) and ([TJJ by 
a linear extrapolation to R^ 1 = of the values obtained 
for at least three different values of R . 

The first and second methods arc compared, for a fixed 
value of the energy E in Fig 12k and for several values of 
E in FigJ^b, and are found to give compatible values of 



the Lyapunov exponent. In practice, the first method has 
the advantage of not requiring a calculation for different 
values of the radius R, the border effects being elimi- 
nated by a suitable choice of the fitting range. Also, the 
value of |t(n)| is bounded from below by the numerical 
accuracy e n um so that the second method may be used 
only if exp(-K-R) > e num . 

A symmetry trick: In practice, to reduce the computa- 
tional effort, we have imposed some symmetry proper- 
ties on the disorder used to obtain the results shown in 
Figs. [T] We have imposed reflection symmetries with 
respect to the planes x = 0,?/ = 0,;z = and we have 
set Tq = 0, so that the independent unknowns in the lin- 
ear system Eq. ([7]) are the amplitudes Di in the sector 
x > 0, y > 0, z > 0. This effectively reduces the number 
of unknowns by a factor of 8, resulting in a gain of a 
factor 64 on the memory size and of a factor 512 in the 
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CPU time, without affecting the vaiue of the Lyapunov 
exponent k, as shown in FigJ3] We have also used this 
symmetry trick in Figs. |4]and[5] 



B. Numerical results for the localization length 

Using the first method described in the previous sub- 
section, we have calculated numerically the Lyapunov ex- 
ponent k and the localization length £ = 1/k for different 
values of the effective scattering length a e s, and of the 
filling factor p occ , see Figs l4|5l6l 

A central point is the choice of the value of p OC c- Our 
scope was to study the regime where the presence of the 
lattice has a small effect on the localization properties, 
since we consider here the lattice mainly as an experi- 
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mental tool to realize strong disorder. To this end we 
will not investigate the regime where p occ only weakly 
differs from unit: This regime indeed corresponds to a 
matter wave propagating in a periodic structure with di- 
lute random vacancies. On the other side, the opposite 
regime p occ — > is not favorable experimentally since it 
leads to large localization lengths scaling at least as the 
mean distance between scatterers d/polc [HI- 



For these reasons, we choose in Fig|4] the reason- 
able value pocc = 0.1. This figure shows the localiza- 
tion length, see Fig0h<, and the Lyapunov exponent, see 
FigHb, as a function of the energy E of the emitting 
source, both for negative and positive values of E. Sev- 
eral values of a ff arc considered in this figure. 
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FIG. 5: (Color online) The same as Fig. [4] except that: The value of p occ is 1/80, the radius is R = 140d (corresponding 
to a mean number of scatterers (TV) ~ 1.4 x 10 5 ), and a e g is rescaled to have the same values of p 1 ^ 3 a e ff as in Fig. [4] where 
p = Pocc/rf 3 is the mean density of scatterers. Since pic? is two times smaller as compared to Fig. [4] we have a e g = 2 x O.lrf 
(black solid line), a e g = 2 x 0.2d (red solid line) , a e g — 2 x 0.7d (green solid line), a e g — 2 x d (blue solid line), a e g = 2 x 1.3d 
(violet solid line). In (a), for 4(E/Eo) = —1.2, one has a cti /d = 2 x 1.3, 2 x 1, 2 x 0.7 from bottom to top, and for A(E/Eo) = 0.2, 
one has a c g/d — 2x 0.7, 2 x 1, 2 x 1.3 from bottom to top. The axes are rescaled accordingly, to allow a direct comparison with 
the results of Fig. [J] The number of realizations of disorder is equal to 100 for all points. 
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The values a c g/d <C 1 are easy to analyze in a per- 
turbativc picture: The matter wave inside the cloud of 
scatterers experiences a mean field shift equal to pg e s 
where p — p occ /d 3 is the mean density of scatterers and 
g c ff = Airti 2 a c g I (2m) is the effective coupling constant 
between the matter wave and a scatterer. In this mean 
field model, the Green's function is simply that of a free 
matter wave in the presence of the uniform mean field 
shift pg c s inside the sphere. We then have, for r < R 
and ro < R 



1 



m e 



ifc m fr 



(16) 



'E + iO+ -Umf ru/ 2ttH 2 r ' 
with "H m f =% + pg c s- The wave vector is then given by 



E 



h 2 k 2 
2m 



h 2 k 2 
2m 



P9cs- 



(17) 



The mean field prediction is thus that, for E < pg e s, k m { 
is purely imaginary, the matter wave cannot propagate 
in the gas of scatterers and is damped with a Lyapunov 
exponent 



K mf = 47rpa off - k 



(18) 



The two conditions of validity for this mean field predic- 
tion are that (i) the matter wave scattering on a single 
scatterer should be in the regime of the Born approx- 
imation, |fca e ff| "C 1, and (ii) there should be a large 
mean number of scatterers within the volume £j^ f , that 
is p£mf ^ 1j so one may neglect fluctuations in the den- 
sity of scatterers over the scale £ m f. The mean field pre- 
diction (TTS)) is plotted as a dashed line in Fig. 2] over the 
range E < pg c s , for the two lowest values of a e ff/d, where 
it is in good agreement with the numerical results. For 
E > pg c f[, the mean field picture cannot of course pre- 
dict how the Green's function decays with |r — ro|. In the 
numerical approach, the localization length turns out to 
be too large to be determined in a reliable way for the 
accessible system sizes. 

In the regime E > p<? G ff , we have found a second man- 
ifestation of the mean field potential pg c g ■ In an optical 
analogy, this potential creates a discontinuity in the re- 
fractive index for the matter wave at the border of the 
sphere containing the scatterers, the index passing from 
unity outside the sphere to n ps (1 — pg c e/ 1 E) 1 / 2 within. 
This induces a reflection of the matter wave, with a Fres- 
nel reflection coefficient r = (n — l)/(n + 1). Since the 
source position ro is in the center of the sphere, this gives 
rise to a radial stationary matter wave with an intensity 
contrast (I max - / min )/(/ mM + I min ) = 2|r|/(l + |r| 2 ), 
provided that one can neglect the attenuation of the co- 
herent part of the matter field, that is if the mean free 
path ps l/(po- scat t), with o- scatt = 47ra^ ff , is larger than 
the sphere radius. We have indeed observed this phe- 
nomenon numerically. 

For larger values of ct e ff, of the order of the lattice spac- 
ing <i, Fig. 2] shows an interesting structure of maxima 



and minima of the Lyapunov exponent at both positive 
and negative energies. 

For E > 0, the value of a e g = 0.7c? shows a wide re- 
gion of remarkably small values of the localization length, 
with £ ps d even smaller than the mean distance between 
scatterers. For larger values of E, a sharp rise of £ is 
observed. A natural question is whether this rise may 
be attributed to the presence of a mobility edge. This 
requires a proof of the presence of localized states to the 
left of this "edge". This will be addressed in section 
TVl where the presence of localized states will be con- 
firmed for a e s/d = 1 and 1.3, in contrast to the cases 
a e s/d = 0.1,0.2 and 0.7. A careful study of the corre- 
sponding phase transition and its critical exponents re- 
quires an examination of finite size scaling that we leave 
for future work. 

Another interesting point is that, further increasing the 
value of a e g/d, although intuitively it may be thought to 
increase the strength of disorder, actually leads to larger 
values of £/d. In particular, at unitarity (a e s/d — > +oo) 
our system sizes were too small to allow a reliable deter- 
mination of £. Similarly, for negative values of a e g/d, we 
have not obtained evidence of a finite £. This is consistent 
with the prediction, based on a pcrturbative calculation 
of the transport mean free path, that no matter wave 
localization can take place for a c s < [l4[ . 

For E < 0, Fig. Q] shows an energy interval where the 
localization length takes very large values. At first sight, 
this may be surprising, as one may naively expect e.g. 
from Eq. (|18p . that k is an increasing function of \E\ 
at least equal to (2m\E\) 1 / 2 /h. In the case where a c fr is 
much smaller than the mean distance between scatterers, 
we relate the existence of this energy interval to the fact 
that the matter wave can form a bound state (a dimer) 
with wavefunction 4>o(r— r^) for each scatterer i, of spatial 
extension a e g and of energy [l9| 



Em 



dim 



2ma 2 cS ' 



(19) 



The various probability amplitudes to have the dimer on 
the sites r^, Tj, . . ., are then coupled by transition ampli- 
tudes itrans of the order of -Edim times overlap integrals 
between </) (r— r.;) and 4>q(y — Tj), integrals that drop ex- 
ponentially with the distance as exp(— |r; — rj|/a ff). A 
precise calculation gives a transition amplitude between 
two scatterers separated by a distance (20| 



>(nj) 



-ry/a ef f 



ma e g 



(20) 



Over the energy interval roughly extending around -Edim 
with a width of the order of typical values of |ttrans|, we 
then face the problem of a bound state that can "tun- 
nel" from one scatterer to another, which may result in 
a localization length £. It turns out that this length may 
have sharp rises as a function of the energy, suggesting 
the occurrence of two mobility edges. This picture can- 
not of course be quantitative for values of a e g as large as 
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0.7 d, but it correctly predicts that the negative energy 
interval of large values of £ shifts to higher energies and 
broadens, when one increases a e s/d. 

In Fig. [5] we consider the same problem as in Fig. 0] 
with a much smaller filling factor p occ = 1/80. We con- 
sidered values of a e g such that p 1 / 3 a e fi is the same for 
the two figures. We also rescaled the axes of the figure 
so that two identical values of the abscissa for Fig. 2] and 
Fig. [5] correspond to an identical value of mE/(p 2 / 3 h 2 ), 
and two identical values of the ordinate correspond to an 
identical value of p 1,/3 £. This allows us to see the effect of 
the underlying lattice [l8j . There are some quantitative 
(not qualitative) differences between the two figures. In 
particular, at negative energies, the separation between 
the energies where £ has a sharp rise is smaller in the case 
of Pocc = 1/80 than in the case of p occ = 1/10 and the 
periodic case p occ — 1 (as we will see in Fig. [6]), on the 
rescaled energy axis. An interesting question is whether 
such a sharp rise is still present in the absence of the 

I 



IV. DENSITY OF LOCALIZED STATES 

Calculation of the localization length, as defined in the 
previous section, is not sufficient to prove the existence 
of localized states at the considered energy. E.g. it is 
expected that, in a spectral gap of the system, the local- 
ization length is finite whereas there is no state available. 
In the present section, we directly investigate the pres- 
ence of localized states as a function of the energy E. 

A. Method 

Wc have to distinguish two cases, according to the sign 
of the energy E. 

Negative energies: For E < 0, the eigenstates of the mat- 
ter wave are discrete bound states. Such a bound state 
corresponds to the matter wave being trapped by the gas 
of scatterers, but it is not necessarily a localized state. 
In particular, in the limit < a e s smaller than the mean 
distance between scatterers, we recall the picture of an 
AB dimer that may tunnel from one scattering site to 
another, see subsection IIII B[ and this dimer, with a neg- 
ative total energy, may be delocalized over the whole gas 
of scatterers. 

The bound state eigenenergies correspond to poles of 
the Green's function Eq.([2|). Since G(r,r ) diverges for 
such an eigenenergy, the secondary source amplitudes Di 
diverge in Eq.©, or equivalently, according to Eq.©, 
this implies that the matrix M defined in Eq.((5]) has a 
zero eigenvalue. For E < 0, the matrix M is real and 
symmetric, so has N real eigenvalues m,. To find the 
bound state energies, there are two strategies. 

In the first strategy, keeping the effective scattering 
length a e & fixed, one numerically calculates the N eigen- 



lattice. 

Finally, we consider in Fig. [5] the "periodic" case cor- 
responding to p occ = 1. This study can be instructive 
to see to what extent the disorder affects the energy de- 
pendence of k. In the periodic case, of course, the Bloch 
theorem forbids the existence of localized states in the 
bulk. The non-zero values of k are then not a signature 
of disordcred-induced localization, but rather result from 
the existence of forbidden energy bands: at forbidden en- 
ergies the matter wave cannot propagate and forms an 
evanescent wave exponentially decreasing away from the 
source. For E < 0, as in the disordered case, an allowed 
energy band is present, which may be interpreted as a 
dimer energy band for a e g much smaller than the lattice 
spacing. For E > 0, the dependence of £ with the energy 
is more flat than in the disordered case, except close to 
the band edges where it diverges. We note that the for- 
bidden band is quite broad for values of a e s of the order 
of the lattice spacing or larger. 

I 

values rrii(E, a e g), 1 < i < N, and then solves for E each 
implicit equation 

mi(E,a eS ) = 0. (21) 

This may be done efficiently by dichotomy, using the 
property that rm is a monotonically increasing function 
of E (see appendix In the second strategy, one in- 
verts the problem, working for a fixed value of the energy 
E and solving mi(E,a e ff) = for a e g. This turns out to 
be straightforward, because of the following structure of 
the matrix M @: 

M = —Id + M 00 (E) (22) 
a c ff 

where depends only on the energy, not on the effec- 
tive scattering length. For each value of E, one then sim- 
ply has to diagonalize (E) , with resulting eigenvalues 
m°°(E). From the relation m,i{E, a cff ) = gQt +m°°(E) 
. one obtains the following parametrization of the neg- 
ative eigenenergy branches: 

— = -mT{E). (23) 

Ocff 

Positive energies: Let us consider first the ideal case of 
scatterers extending over the entire position space. Then, 
for E > two kinds of states are expected, square inte- 
grable localized states corresponding to a discrete spec- 
trum, and extended states corresponding to a continuous 
spectrum. It is expected that the continuous spectrum 
exists at E > E c , where E c is called a mobility edge, and 
that the localized states are essentially at E < E c |22| . 

In reality, the scatterers occupy a finite volume of ra- 
dius R, so that, for E > 0, the energy spectrum of the 
matter wave forms a continuum extending from to +oo, 
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FIG. 6: (Color online) The same as Fig. 21 except that: The value of p occ is 1 (there is no disorder, since all the lattice sites 
within the sphere of radius R = 30d are occupied by the scatterers, which leads to N ~ 1.1 x 10 5 [IJ), and a e g is rescaled 
to have the same values of p 1 ^ 3 a e s as in Fig. [4] where p — p OC c/d 3 is the mean density of scatterers. Since Pocc is (O.l) 1 ^ 3 
times smaller as compared to Fig. [4] we have o e g = (0.1) 1/3 x O.ld (black solid line), a cB = (0.1) 1/3 x 0.2d (red solid line) , 
a cS = (0.1) 1/3 x 0.7d (green solid line), a e s = (0.1) 1/3 x d (blue solid line), a e g = (0.1) 1/3 x 1.3d (violet solid line). In (a), for 
{0.1) 2/3 E/E Q < -1.2, one has a cH /d = (0.1) 1/3 x 1.3, (0.1) 1/3 x 1, (0.1) 1/3 x 0.7 from bottom to top, and for (0.1) 2/s E/E = 0.5, 
one has a cB /d = (0.1) 1/3 x 0.7, (0.1) 1/3 x 1, (0.1) 1/3 x 1.3 from bottom to top. The axes are rescaled accordingly, to allow a 
direct comparison with the results of Fig. 2] 



and none of the positive energy eigenstates is square in- 
tcgrable. Nevertheless, for a large enough radius R, one 
can still find eigenstates which, inside the medium, expo- 
nentially decrease over several orders of magnitude away 
from a central region, toward the borders of the gas of 
scatterers @. For all practical purposes, these are local- 
ized states. From a spectral point of view, these localized 
states are expected to correspond to very narrow reso- 
nances of the matter wave inside the scattering medium, 
that is to complex poles 

z rcs = E TCS - ihT/2 (24) 

of the Green's function Eq.([2]) analytically continued to 
the lower half of the complex plane, T > 0. Such a res- 
onance state, localized close to the center of the sphere 
containing the scatterers, has a width T, i.e. an inverse 
lifetime, that vanishes exponentially with the radius R 
as exp(— 2R/£), where £ is the localization length of the 



localized states 0, [l4j], reproducing for R = +oo the 
aforementioned ideal case of poles on the real axis. 

To find such resonances, one has to analytically extend 
the Green's function Eq.Q to complex energies z. Then, 
setting k = (2mz/h 2 ) 1 ^ 2 in Eq.©, taking the real and 
negative axis as the line cut for z 1 / 2 , one finds that Eq.© 
still holds. Then the complex poles z les are such that the 
matrix M has a zero eigenvalue. Adopting the second 
strategy used for bound states, one actually diagonalizcs 
M°°(z) so that one has to solve the equations for z rcs : 

— = -mHzres). (25) 

O e ff 

In practice, to solve Eg. (f2"5]l , we implement Newton's 
method as follows. We choose some real and positive 
value for E. We diagonalize M°°(E), calculating the 
eigenvalues and the eigenvectors; we select an eigenvalue 
m°°(E) and we set the effective scattering length to the 
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value such that a 



\E) is purely imaginary, 



1 



a c ff = 



Rem°°(E)' 



(26) 



The successive steps of the method involve the calcula- 
tion of the derivative of m°°(z) with respect to z, which 
can be done thanks to the extension of the Hellmann- 
Feynman theorem to non-hermitian matrices [23j . 



This constitutes the initial guess in Newton's method. 



J 



Of particular interest is the approximate value of z p 
obtained after the first step: 



E 



Im mf(E) 
' dm°°(E)/dE : 



(27) 



which, if compared with (|24[) . provides aproximate values 
of E Tes and T. We have checked that, if the imaginary 
part of m°°(E) is small enough, typically < 10 -3 /<i for 
Pocc = 0.1, the first step Eq.(f2"T|) has in practice already 
converged to the exact pole location. To obtain FiglJl we 
used this first step approximation. 



B. Results for the density of bound states and 
resonances 

Using the techniques presented in the previous subsec- 
tion, we have calculated the density of resonances and 
of bound states in the plane (E,a e g), E being the real 
part of the Green's function complex poles for resonances 
and the eigenenergy for bound states. The raw result is 
presented in Fig(7k.. 

For large and negative E, and positive a c g, a first class 
of bound states is observed, with a density of states con- 
centrated in a narrow interval of values of a c g. This 
was expected: For a e g > the matter wave can form a 
bound state (a dimer) with a single scattercr, of energy 
Edim given by Eq. ([T5|) . In the presence of the disordered 
ensemble of scatterers, this dimer subsists in the limit 
Ocff "C p -1 / 3 , where p is the mean density of scatterers. 
Simply, the dimer state may tunnel from one scattercr 
in location to another one in Yj , with a tunneling am- 
plitude close to Eq. ([2Hl) . This tunneling broadens the 
eigenenergy interval around .Edim by an amount of the 
order of the tunneling amplitude, a small amount in rel- 
ative value if a c s ^ p^ 1 ^ 3 . 

As one moves to larger values of the energy, the afore- 
mentioned interval of eigenenergies moves toward larger 
values of a e g and broadens, as expected from the pre- 
vious reasoning. In addition, an internal structure ap- 
pears in the figure, see the zoom in FigJHl the density of 
states being peaked over several subintervals, correspond- 
ing to various branches. These branches may be identi- 
fied physically by calculating the eigenenergies of trimers, 
tetramers, ... of the matter wave with two, three, ... scat- 
terers on the lattice, see the lines in the corresponding 
FigJHl One simply has to apply the method described 
in the previous subsection, for N = 2, 3, . . ., sec Eq. 



and the corresponding matrices M(E) may be diagonal- 
ized analytically. Setting E = —h 2 q 2 /(2m),q > 0, one 
finds that a e g is given as a function of the energy of an 
AB-2 trimer by the two branches 



1 e~ qri2 

=q± 

a e g r 12 



(28) 



where r\2 is the distance between the two B scatterers. 
The calculation may be done also for the AB3 tetramers, 
giving rise to three branches for each particular set of val- 
ues of the inter-scatterer distances ria, T13, ^3. The an- 
alytical expressions are simple in the relevant case 7*13 = 
r 2 3, where one sets /3 = e~ gri2 jryi and 5 = e~ 9ri3 /ri 3 . 
One then finds that one of the branches reduces to one 
of the trimer branches, 1/gmj = Q + /3. The other two 
branches are then given by fl9l ]: 



1 

floff 



<1 



/?±(/3 2 +85 2 ) 1/2 



(29) 



The values of a e g for trimer and tetramer thresholds are 
readily obtained from (|28| and (|29|) by taking q — > 0. 
Remarkably some of these threshold values are negative. 

For negative values of E close to the origin, another 
class of bound states is observed for negative a c g. For 
values of a e g close to zero, this class may be interpreted 
by a simple mean field effect. The matter wave experi- 
ences an effective attraction inside the volume containing 
the scatterers, represented by the negative mean field po- 
tential pg c fi , see Eq. (fTT|) . so that the gas of scatterers pro- 
duces an effective square well potential of radius R that 
can support matter wave bound states. This mean field 
interpretation is in good agreement with the numerical 
results for small | a c ff | . The oblique mean field white line 
in Fig(7] giving the bottom pg c s of the effective square 
well potential accurately reproduces the lower boundary 
of the bound state energies. We note that expressions 
(|2"8")l and (|2"9"1) predict the presence of bound states also 
at negative values of a e g, for low enough values of q. In 
such an energy region the values of a e g are large and 
negative (a c g < —d/2). Then the effective scattering 
length is typically larger that the mean distance between 
scatterers p -1 / 3 < |a e ff|- This may explain why no well 
delineated internal structure is observed in the density of 
states for negative a e g. 

Each of the two components that we discussed for 
E < continuously develops in two corresponding com- 
ponents on the E > side of Fig[7^,. 
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(a) (b) 




FIG. 7: (Color online) For the 3D system: Density of resonances and bound states per scatterer in the plane (energy E, 
effective scattering length a e g), obtained as explained in subsection IIV Al E on the horizontal axis is either the real part of 
z ICS on the positive energy side, or the bound state eigenenergy on the negative energy side. The filling factor is p occ = 1/10 
within a sphere of radius R = 29d, so that the mean number of scatterers is (A) ~ 10 4 . The value of E is discretized with 
a step 0.0025i5o for both positive and negative values of E, and we use Eo = ti 2 /(md 2 ) as unit of energy. For each value 
of E one different realization of disorder is used, without imposing any reflection symmetry (see end of subsection IIII A|) . A 
logarithmic scale is used: The color map (see bar on the right) is applied to the quantities log 10 — j^p-^g , for E > 0, and 
log 10 i^apd.^, for E < 0, where N leB and A bound are the number of resonances and bound states respectively, within each 
rectangular bin of area 8S = 8ESa g (SE = O.OliJo and da e g = 0.007d, for both positive and negative values of E). The plotted 
quantity is then ~ log 10 (1.4A ros ), and ~ log 10 ( 1.4 Abound), respectively. The horizontal dashed lines correspond to the values 
of a e g/d £ {0.1,0.2,0.7, 1, 1.3} used in Fig. 3] The oblique solid line is the border of the energy gap E = pg e ff predicted by 
the mean field theory (see Eq. (|18|l with « m f = 0). (a) No selection is applied to the resonances (E > 0). As explained in the 
text, most of the displayed resonances (the ones with too large a width) are not expected to be meaningful, (b) Restricting to 
E > 0, only the resonances with a width T < r max = 10~ 6 Eo/h are kept in the density of resonances. The value of F ma x is 
essentially infinite with respect to the duration of typical experiments: For a matter wave of 87 Rb atoms and an optical lattice 
with d = OA/jm, as in [ll[, one has indeed l/r max ~ 200 seconds, (c) Only the resonances and bound states corresponding to 
a small enough participation volume V p are kept in the density [Vp 1 /d < 6.5, see Fig(9^ and c]. (d) Only the resonances and 
bound states corresponding to a small enough r.m.s. size a in real space are kept in the density (cr/d < 4.2, see Fig|5p and d). 
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FIG. 8: (Color online) Zoom of figure [7^., in the region of 
E < and small and positive values of a c a (here 8E = 0.01 Bo, 
and 8a e B = 0.005(f). On the color map representing the den- 
sity of states (see caption of Fig[7£i), are plotted the ener- 
gies of bound states of an A atom with one, two or three 
B atoms, denoted respectively dimer AB, trimer AB2, and 
tetramer AB3 states. More precisely are plotted the energies 
of a dimer Eq. (|19[) (black solid line), of a trimer Eq. (|28|) with 
two B atoms separated by a distance ria = d (dashed black 
and dash-dotted black lines), of a trimer with two B atoms 
separated by a distance na = V%d (dashed and dash-dotted 
blue lines), and of a tetramer Eq. (|29[) with three B atoms sep- 
arated by distances ri3 = r23 = d and na = y/2d (solid and 
dashed magenta lines). It is worth noting that the prediction 
for the energy of such isolated bound states correspond to the 
higher density lines of the color map plot. 



in Fig[7j3. the component in the lower part of FigJ?^, con- 
tains no long lived resonances, hence no localized states. 

The E > component in the upper part of FigJTJi is 
more promising for observing localized states. In FigfTjD 
we have kept only very long-lived resonances, with a de- 
cay rate T < 10~ 6 h/(md 2 ) [24j]. This corresponds to a 
lifetime larger than ~ 200s for 87 Rb atoms in an opti- 
cal lattice of spacing d = 0.4^m [ll[, and can then be 
considered as infinite with respect to the typical dura- 
tion of ultracold gases experiments. We see that, after 
this selection, the upper component still gives a signifi- 
cant contribution, while the lower component essentially 
disappears. 

In FiglZfc and Fig|7pl, we have kept only the bound 
states and resonances that are localized in real space, 
according to the two following criteria respectively: In 
FigEb we keep states with a participation volume V p 
smaller than (6.5c?) 3 , where the participation volume is 
defined as: 



(30) 



where the amplitudes Di are the components of the eigen- 
states of the M matrix with a zero eigenvalue, and are 
normalized such that Ej=i |A| 2 = 1- The participation 
volume (also called participation ratio) provides an esti- 
mate of the volume over which the amplitudes Di take 
significant values In Fig|7ji we keep states with a 
root-mean square size a in real space smaller than 4. 2d, 
with the definition 



J2 — 



< N 



N 



I A 



(31) 



The E > component in the lower part of Figj7^ has 
a first boundary (delimiting a blue triangular zone con- 
taining no resonances) which is well described by mean 
field theory at low values of E and a c ff, see the oblique 
white line. It has a second boundary on the right, corre- 
sponding to a blue oval region again with no resonances. 
There is no physical interpretation of this blue oval region 
at the present moment. A word of caution may be useful 
at this stage. For E > 0, we calculate a density of reso- 
nances, with an approximate method valid for long-lived 
resonances only (see subsection IIV Ap . so the absence of 
resonances does not imply the absence of eigenstates (in 
particular extended states) of the system. The blue oval 
region thus does not necessarily correspond to a spectral 
gap in the density of states. In particular, the calculation 
of the localization length £ performed in the previous sec- 
tion for a e e/d = 0.1 shows a divergence of £ at the mean 
field border, with small values of £ at energies below this 
border, which is compatible with a spectral gap. However 
it does not show any sharp decrease of £ at the border 
with the oval region, £ retains very high values even for 
energies in the oval region, raising doubts about the ex- 
istence of a spectral gap in that region. As we shall see 



where the Di are defined as for V p . 

For computational convenience, in defining these two 
quantities we used the amplitudes Di, rather than the 
wavefunctions tp of the resonances and the bound states. 
We verified for several examples that localization of Di 
implies a localization of tp, and vice versa. This can be 
easily understood since the matter wave problem is anal- 
ogous to that of a scalar wave of light scattered by pinned 
atomic dipoles. In this case, the values of Di correspond 
to the values of the electric dipole moment on each atom, 
and the matter wave wavefunction ijj corresponds to the 
electric field £. In this scalar light analogy it is intuitive 
that the spatial extension of the dipoles reflect that of 
the electric field, and vice versa [T^| . For completeness, 
we provide the expression of ip in terms of Di in the ap- 
pendix [B] 

The choice of above mentioned limiting values for V p 
and a are motivated by the histograms of FiglHl For 
E > 0, the histograms reveal a bimodal structure (i.e. 
two maxima) in the probability distribution of V p and 
a (the bimodal structure is quite clear in Fig. [SJi, and 
less clear in Fig. [3b), which suggests the coexistence 

1 /3 

of localized resonances (small values of V p and a) 
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FIG. 9: For the 3D system: For all resonances with E Tea /Eo € 
(0, 2) and a e g/d G (—3, 3), and all bound states with E/Eq € 
(—2,0) and a a g/d £ (—3,3), the figures show an histogram 
giving the number of resonances (in (a) and (b)) and bound 
states (in (c) and (d)) per scatterer (the number of scatterers 
is N = 10 4 ) as a function of: in (a) and (c), the cubic root of 
the participation volume defined in Eq. (|30[) ; in (b) and (d), 
the r.m.s. size in real space defined in Eq. (|31(l . The bin size 
is ~ 0.2d for (a) and (c), and ~ 0.26d for (b) and (d). The 
parameters are the same as in Fig. [7^,. The dashed vertical 
lines are the values Vp 173 /d ~ 6.5 and a/d ~ 4.2 [2^| used in 
Fig. [7Jj and Fig. [TJi to select the bound states and resonances 
that are spatially localized and to filter out the bound states 
and resonances that are spatially extended. 



and extended resonances (large values of V p and a, 
of the order of the radius R = 29d of the sphere con- 
taining the scatterers). For E < 0, the bimodal struc- 
ture still appears for a, showing that some bound states 
are extended. Interestingly, the bimodal structure is not 

1 /3 

present for E < on the histogram oiVp'. Actually, it 
may happen that a state which is a coherent superposi- 
tion of a few-body bound state (e.g. a dimer, a trimer) 
at different locations is considered as a localized state 
with the criterion based on V p , whereas it is considered 
as an extended state with the one based on a. An ex- 
treme example of such a case corresponds to a dimcr 
delocalized over two sites of positions r and — r. Then 

Di = (<5 ri , r + S rz ^ r )/V2, leading to V^ 73 = (2/p) 1 / 3 
(insensitive to the distance 2r between the two possible 
dimer locations), and to a = r [26j . 

Coming back to FigfTfi, wc thus see two main streams 
of localized states, one restricted to negative energies, 
and the other extending to positive energies. The exis- 
tence of these two streams explains the structure of the 
localization length in FigfJJt. It is possible to establish a 
correspondence between the occurrence of divergences of 
the localization length, and boundaries between an en- 
ergy interval with localized states and an energy interval 
with extended states. A particularly rich example corre- 
sponds to the effective scattering length value a e g = d, 



for which by combining the information coming from the 
localization length £ in FigfJJt, and the density of reso- 
nances/bound states in FigsfTk and [7H, we can distin- 
guish 5 energy intervals: 

• E/Eq < —1.3: the density of bound states is zero 
in FigfT^, and £ takes small values in FigJJJt. This 
reflects the occurrence of an energy gap. 

• —1.3 < E/Eq < —1: the density of localized bound 
states is non-zero in FigJTJl, and £ increases with a 
divergent behavior at E/Eq ~ — 1 in FigfJJt. This 
suggests the presence of an energy band of localized 
states with an upper mobility edge at E/Eq ~ — 1. 

• — 1 < E/Eq < —0.2: there are bound states in 
FigEK but no localized bound states in FigfTH, 
and £ assumes very large values (presumably larger 
than the system size). This suggests the presence 
of an energy band of extended bound states de- 
limited by two mobility edges at E/Eq ~ —1 and 
E/Eq ~ -0.2. 

• —0.2 < E/Eq < 0.6: the density of localized 
bound states and localized resonances is non-zero 
in FigfTJi, and £ assume small values away from di- 
vergent behaviors at E/Eq ~ —0.2 and E/Eq ~ 0.6 
in FigfJJt. This suggests the presence of an energy 
band of localized states delimited by two mobility 
edges at E/Eq ~ -0.2 and E/Eq ~ 0.6. 

• E/Eq > 0.6: the density of localized resonances is 
zero in FigJTJi, and £ assumes very large values (pre- 
sumably larger than the system size) . This suggests 
the presence of an energy band of extended states 
delimited by a lower mobility edges at E/Eq ~ 0.6. 

One can apply the same analysis to other values of 
a c fr considered in FigJJJt. In particular, it is interesting 
to note that for a e ff = 0-7d the correspondence between 
the positions of mobility edges deduced from the local- 
ization length £, and deduced from the density of res- 
onances/bound states, is only qualitative: for negative 
energies, £ diverges at E/Eq ~ —1.2 and E/Eq ~ —0.8 
in FigfJJt, while one finds the transition region between 
localized and extended states at E/Eq ~ —1.0 and 
E/Eq ~ —0.9 in FigfTfl. This discrepancy can be ex- 
plained in terms of finite size effects, as appears from 
FiglHH Indeed, FigfTUk shows a zoom of the figure [7H 
(i.e. density of resonances and bound states per scatterer 
in the plane (i?,a e ff) having a r.m.s. size in real space 
a < 4. 2d). In FigllOb wc impose an additional filter- 
ing with respect to that of FigfTUk: we keep only bound 
states and resonances having their maximum |_D,-| at a 
distance less than R' = 20d from the origin, whereas the 
disorder fills a sphere of radius R = 30d. In this way only 
localized states located far from the border of the system 
are selected. In figure [TOb the transition region between 
localized and extended states is now at E/Eq ~ —1.2 
and E/Eq ~ —0.8, in quantitative agreement with the 
mobility edges deduced from £ in FigfJJt. 
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FIG. 10: (Color online) For the 3D system: Finite size effect 
on the density of localized states, (a) Zoom of the figure [7ji, 
which shows the density of resonances and bound states per 
scatterer in the plane (_E,a e ff) having a r.m.s. size in real space 
a < 4.2rf. We recall that there are N ~ 10 4 scatterers in a 
sphere of radius R = 29d. (b) Same as (a), with an additional 
filtering: we keep only bound states and resonances having 
their maximum \Di\ at a distance less than R' = 20d from 
the origin (in the normalization, we considered the number of 
atoms N' ~ 3.3 x 10 3 in the sphere of radius R'). 



useful for experimental purposes to identify optimal val- 
ues of {E, cieff) for the observation of maximally localized 
resonances. These occur for values of a g 0.8d and 
< E/Eq < 0.4, for the considered density of scatterers 
Pocc = 1/10. 

V. LOCALIZATION IN A 2D GEOMETRY 

It is now experimentally possible to realize systems 
where the atomic matter wave species A is confined in a 
2D geometry. One can apply a strong confinement along 
the direction z, freezing the motion in its ground state 
along z, while leaving a free motion in the xy plane. In 
analogy with the 3D disordered system studied in previ- 
ous sections, we study here the 2D case, where both B 
scatterers and matter wave atoms A are confined along 
the z direction (not necessarily with the same potential) , 
and B scatterers are randomly distributed at the nodes 
of a 2D optical lattice in the xy plane (the A atoms being 
insensitive to that lattice). As for the 3D case, the low 
energy scattering of an A atom with a single trapped B 
atom can be characterized by a 2D effective scattering 
length, accounting for the effect of the atomic motion of 
the B atom during the scattering process. This allows us 
to replace the B atoms by point-like scatterers at fixed 
positions, described by contact conditions on the A atom 
wavcfunction. For simplicity we use for the effective scat- 
tering length in 2D the same notation a e ff as in 3D. To 
our knowledge, the dependence of the effective 2D scat- 
tering length on the free space 3D scattering length a, on 
the atomic masses and on the oscillation frequencies, has 
not yet been investigated. 

On the theoretical side there has been significant in- 
vestigation of disordered systems in 2D since this is the 
lower critical dimension for the occurrence of the metal- 
insulator transition. In 2D it is expected that all states 
are localized in an infinite system, but the localization 
length increases rapidly with the energy 0, [27|. It is 
thus important to identify the range of parameters (mat- 
ter wave energy E and a e g) for which the localization 
length £ takes small values. In subsection IV Al we explain 
how to adapt the 3D formalism to the 2D case, and in 
subsection I V Bl we present the results for the localization 
length and the density of localized states. 

A. Formalism for 2D systems 



To conclude the analysis of the 3D system, we give in 
Fig fill the values of the width r for the resonances in 
the plane (E,a e g). A filtering was applied to these reso- 
nances: In Fig lllb , we considered only resonances with a 
participation volume V p jd < 6.5, and in FigfTTb only 
resonances with a r.m.s size in real space a jd < 4.2. This 
analysis shows the presence of numerous extremely long- 
lived resonances with T < 10 -6 Eq / H. This figure may be 



In 2D the contact condition on the matter wave wave- 
function in the vicinity of each scatterer of position 
is 

V^ r ) = -£j AM|r-ri|/a eff ) + 0(|r- ri |) (32) 
7m z 

where a e g is now the 2D effective scattering length (al- 
ways positive) [28]. Comparing to the 3D contact con- 
ditions Eq. ([1} , we see that — 1/ ln(a e fi / d) in 2D plays 
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FIG. 11: (Color online) For the 3D system: Width F of the 
resonances as a function of the energy E and the effective 
scattering length a e a [see Eqs. (1241) and (|27[) ]. The physical 
parameters are the same as in Fig. [7] The plane (E, a e ff) is 
decomposed in rectangular bins of widths SE and 8a c g. The 
color map (see bar on the right) is applied to the quantity 
log 10 h{F) /Eo where (r) is the mean value of V for the reso- 
nances within a given bin. The resonances are filtered in (a) 
over the participation volume as in Fig. 0: and in (b) over 
the r.m.s. size a as in Fig. 0i. In (a) one has 8E = 0.009-Eo 
and Sa c ff = 0.02d, and in (b) one has SE — 0.008i?o and 
5a eB = 0.02d. 



the role of a e g/d in 3D. We recall that, in 2D, the limit 
Q-off does in fact correspond to a weakly repulsive 
limit, while a e g — > +oo corresponds to a weakly attrac- 
tive limit. 

The 2D Green's function still obeys Eq.Q, where A r 
is now the 2D Laplace operator, so that formally Eq.([6j 
still holds provided that one takes for go the 2D free 
matter wave Green's function. For E > 0, we set E = 



H k ■ /(2m). k > 0, and we obtain 



(33) 



with the Hankel function expressed in terms of Bcsscl 
functions as Hq 1 \z) = Jq(z) + iN$(z). From [29| one 
obtains the limiting behaviors 



H^(kr) 
r(i 



r— »-(-oo 



rfcr 



1/2 



i(kr-ir/4) 



(34) 



K'(kr) = l + -ln(^)+ (l), (35) 

where 7 = 0.57721566 ... is Euler's constant. The nu- 
merical factor in Eq. ([33)) thus results from the fact that 
A r lnr = 27T(5(r) in 2D. For E < 0, we set E = 
— fi, 2 q 2 /(2m), q > 0. That is, k = iq is now purely imagi- 
nary, and we obtain 



5o ( r ) = M K o(<l r ) 



(36) 



where K$ is a modified Bessel function of the second kind. 

The secondary source amplitudes Di are then deter- 
mined by imposing the contact conditions (|32[) at the 
order 0(1). That is, for the non-diverging term lna e ff, 

N , 2 

^ MijDj = —^—g {vi — r ), Vi G {1, . . . , N} (37) 



where we have introduced the N x N matrix 



Ma- 



goiTi-Tj) if lytj, 

m 



~i? +ln(^) if i = j. 



(38) 



For E < 0, where k = iq, q > 0, this holds with the de- 
termination In i = iir/2, leading to the diagonal element 

For the calculation of the localization length £, we use 
the second method already discussed for the 3D case. By 
following the same derivation as in 3D, and using Eq.tfM]). 
we obtain the same expression (fT2"j) for the values of the 
transmission amplitude t(n) of a source located at the 
position r . Then, we define the Lyapunov exponent n 
as 



2tt 



d6 (ln|t(n)|) 



lim 

r-^+oo J 2tt R 
lim / 



for E > 0, (39) 



dd (ln|i(n)|) 



2vr 



for E < 0, (40) 



where n has coordinates (cos 9, sin 9) in the xy plane, and 
R is the radius of the disk containing the disorder. 
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To obtain the density of bound states and resonances comes from setting a e g = d in Eq. ([38|) . To evaluate the 



in 2D, we use the same procedure as in 3D, and split the 
matrix M as 



M 



ln(a e ff/d)Id 

where M x is the value of M for 1 / ln(a e ff /d) 



(41) 



density of bound states for a fixed value of a e g one can 
solve by dichotomy the 2D equivalent of Eq. ([21]) . using 
the fact that the eigenvalues of M(E), for E < 0, are 
a monotonically decreasing function of E (see appendix 

E}. 





E/E Q 

FIG. 12: For the 2D system: localization length £ (a), and Lyapunov exponent k (b), as a function of the energy E. Here the 
2D effective scattering length is a c g = d, the lattice filling factor is p occ = 1/10, and an average over 100 realizations for the 
positions of the scatterers has been performed. The value of k is obtained using the second method [see Eqs. (|39[) and (I40[) ]. 
For < E/Eo < 0.24, we considered larger and larger system sizes 2R/d £ {100, 150,250,350} until a convergent value for k 
is reached, within the statistical error bars, which are smaller than 5%. For 0.24 < E/Eq < 0.5 we considered system sizes 
2R/d £ {350, 600}, and found a maximal deviation of 10% in k between the two sizes, finally plotting the values for 2R/d — 600. 
For negative energies —2 < E/Eo < we considered system sizes 2R/d £ {100,250} (except for E/Eo = —0.02 and —0.01, 
where we considered 2R/d £ {50, 75, 100}), from which we extrapolated linearly to the case 1/R —> 0. In this region we obtain 
error bars smaller than 1%. The absence of the E/Eo = point is due to the fact that the matrix M is logarithmically divergent 
for such a value. 



B. Results for the localization length and the 
density of states for the 2D system 

In figures [T2"a and b wc show the value of the local- 
ization length £ and of the Lyapunov exponent k = l/£ 



of Eqs. ([39j) and (j40|) . as a function of the energy E of 
the matter wave, for the 2D effective scattering length 
«cff = d. For E > we note a monotonically increasing 
behavior of £, and a monotonically decreasing behavior 
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FIG. 13: For the 2D system: Lvapunov exponent k of Fig |12b . 
plotted in logarithmic scale (solid black line). The prediction 
of the independent scattering approximation (ISA) theory is 
also shown (dotted black line). 



of k with the energy E, with the occurrence of a flatter 
region (a "shoulder") around E/E = 0.05. The logarith- 
mic plot for k in FigUjJ] suggests an exponential decay of 
k at large energies. However, being over less than one 
decade, it is not fully conclusive. We also performed cal- 
culations of k for a c g/d = 0.3, obtaining similarly large 
values of k at low energy. In contrast, for a e s/d = 3 the 
values of k drop by an order of magnitude for the same 
energy interval. In FigQ~3] we also plot the prediction 
coming from the frequently used independent scattering 
approximation (ISA) (see, for instance [30[), according 



to which the localization length £ = lexp (irl~Rck c ff / 2) , is 
expressed in terms of the mean free path 1=1/ (2Imfc e gf), 
and an effective wavevector k C R defined as 



h 2 k 2 



^ > p(k\T(E + l + )\k) 



(42) 



2m 2m 

where T is the T-matrix in 2D for a single scatterer [3l| 



(k|T(£ + z0+)|k) = 



-Tvh 2 /m 



ln(fca cff eT/2) - in/2' 



(43) 



The figure shows that the ISA is not accurate for the 
considered parameters, i.e. for large values of a s, of the 
order of the mean distance between scatterers. 

In figure [14] we show the density of bound states and 
resonances in the (E/Eo,]n(a e s/d)) plane, in a 2D sys- 
tem (see sections llV Al and IV A[) . Calculations have been 
performed for a disk of radius R = 150d, and a den- 
sity of scatterers p occ = 1/10. In FigfMk we show raw 
data without any selection, in Fig |14b we present only 
resonances with a width V < 10~ 6 Eo/h, in FiglMb we 
present only resonances and bound states with a partic- 

1 /2 

ipation surface S p /d < 9.5 (this value is motivated by 
FigsfTTk and c, see below), and in FigfT4"H we present only 
resonances and bound states with r.m.s. size a/d < 4.2 
(this value is motivated by Figsfl7b and d, sec below). 
Here the participation surface is 



(44) 



where the Di are normalized as y \D/\ 2 = 1, and a is 
defined by the 2D version of Eq. (|3"Tj) . 



In figure [TBI we show a zoom of FiglMk in the region of 
lowest energies. The streams of higher densities of bound 
states are shown to be in correspondence with the energy 
of some few body bound states (AB dimers, AB2 trimers, 
AB3 tetramers). We recall that, in 2D, the matter wave 
has a single bound state AB on a isolated scatterer, with 
an energy pjl, [H| 



•Edim — — 



2h 2 



ma 2 cS e 2 T ' 



(45) 



The energy Strim 
given by [19( 



h q /(2m) of the AB2 trimer is 



d c ff 1 

— - — = - cxp [±K (qri 2 )], 
A q 



(46) 



where r\% is the distance between the two B scatterers. 
The energy EWra = —h 2 q 2 /(2m) of the AB 3 tetramer 
is given [19j . in the particular case of ri3 = r23, by the 



expression 



a off e 



7 



■ exp 



-W±((3 2 + 5 2 y/ 2 ] 



(47) 



where we only considered the branches -Etetra 7^ E tr i m , 
(3 = K (qr 12), 5 = K (qr 13 ) = K (qr 2 3), and r„- is the 
distance between the two B scatterers i and j. The signa- 
ture left by few body bound states in the density of states 
can be further appreciated in Fig[16l where we plot the 
density of states as a function of the energy E < 0, at the 
fixed value of a c ff = d. The energies of some few body 
bound states (vertical arrows in the figure) are shown to 
correspond to peaks in the density of states. Other peaks 
can be identified with other few body bound states (not 
shown). The same calculation, performed for a G ff = 3d, 
shows a broad structure where essentially all the peaks 
are washed out. 

In figure fTTl we show histograms of the number of states 
as a function of the participation surface S p and of the 
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FIG. 14: (Color online) For the 2D system: Density of resonances and bound states per scatterer in the plane (energy E/Eq, 
logarithm of the effective scattering length ln(o e fj/d)), obtained as explained in the subsection IV Al E on the horizontal axis is 
either the real part of z ICS on the positive energy side, or the bound state eigenenergy on the negative energy side. The filling 



factor is p OC c = 1/10 within a disk of radius R = 150d, so that the mean number of scatterers is (JV) 
of E is discretized with a step 0.0025-Eo for both positive and negative values of E, and we used Eq - 



; 7 x 10'. The value 
h 2 /(rnd 2 ) as unit of 

energy. For each value of E one different realization of disorder is used, without imposing any reflection symmetry (see end of 
subsection MI A[) . A logarithmic scale is used: The color map (see bar on the right) is applied to the quantities log 10 ^j^-j^fi 



E d 



for E > 0, and log 10 — - N jg- 

within each rectangular bin of area 5S 



for E < 0, where N TeB and Abound are the number of resonances and bound states respectively, 
SES\n(a eS /d) (SE = 0.01E and 5\n(a cS /d) = 0.007, for E > 0, and SE = 0.01£ 



and Sln(a c fi/d) = 0.004 for E < 0). (a) No selection is applied to the resonances (E > 0). As explained in the text, most of the 
displayed resonances (the ones with a too large width) are not expected to be meaningful, (b) Restricting to E > 0, only the 
resonances with a width F < r max = 10~ 6 Eo/h are kept in the density of resonances. The value of F max is essentially infinite 
with respect to the duration of typical experiments, (c) Only the resonances and bound states corresponding to a small enough 
participation surface 5^ are kept in the density (Sp^ 2 /d < 9.5, see Fig |17b . and c). (d) Only the resonances and bound states 
corresponding to a small enough r.m.s. size a in real space are kept in the density (cr/d < 4.2, see Fig |17b and d). 
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r.m.s. size a, for positive energies (FigsflTk and b) and 
for negative energies (FigslTTb and d). The structure of 
the histograms at positive energies shows broad struc- 
tures corresponding to states extending over the whole 
disordered system. This fact does not contradict the gen- 
eral statement that all states are localized in a 2D system 
[27| . Indeed, at large energies and large values of a e e, 
the localization length £ becomes very large, larger than 
the system size. For negative energies, the histograms 
show one narrow peak with no tails, indicating that all 
the states are localized with a localization length of few 
lattice spacings only. 
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VI. CONCLUSION 

We performed a quantitative study of the 3D and 2D 
strong localization of matter waves in a random poten- 
tial realized by point-like scatterers (atoms) pinned at 
the nodes of a cubic or square lattice. This model al- 
lows for an exact numerical analysis of both the local- 
ization length and density of states, as functions of the 
matter wave energy and of the effective scattering length 
between the matter wave and a single scatterer. We con- 
sidered systems having a number of scatterers of the same 
order as that achieved in current experiments with ultra- 
cold atomic gases (N ~ 10 5 ), corresponding to systems 
with a diameter ~ 140 lattice spacings in 3D for a lattice 
filling factor = 1/10. 
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FIG. 15: (Color online) Zoom of figure 114b . in the region 
of E < (here SE = 0.01E , and S\n(a eS /d) = 0.003). On 
the color map representing the density of states (see caption 
of Fig |14b h are plotted the energies of bound states of an 
A atom with one, two or three B atoms, called respectively 
dimer AB, trimer AB2, and tetramer AB3. In particular are 
plotted the energies of a dimer Eq. (|45[) (black solid line), of 
a trimer Eq. (|46p with two B atoms separated by a distance 
f 12 = d (dashed black and dash-dotted black lines) , of a trimer 
with two B atoms separated by a distance r\2 = y2d (dashed 
and dash-dotted blue lines), of a tetramer Eq. (|47p with three 
B atoms separated by distances ri3 = r-23 = d and ri2 = V2d 
(solid and dashed magenta lines). It is worth noting that 
the predictions for the energies of such isolated bound states 
correspond to the higher density lines of the color map plot. 



In figures fT8k and b, we plot the value of the width 
r of the resonances, in the (E/Eo,ln(a e s/d)) plane. In 
FigfTBh,, we considered only resonances with a partici- 

1 /2 

pation surface S p jd < 9.5, while in Fig|18b. we con- 
sidered only resonances with a r.m.s size a/d < 4.2. In 
both figures wide regions are present corresponding to 
very small values of T, i.e. to extremely long-lived res- 
onances (in practice infinitely long-lived at the scale of 
the experiments). 



FIG. 16: (Color online) For the 2D system: density of bound 
states for a given value of the 2D effective scattering length 
deff = d. The histogram is the result of the average over 100 
random realizations of the scatterers positions in a disk of 
radius R = 50d, and a filling factor p occ = 1/10. On average, 
we find 0.85 bound states per scatterer. As done in Fig |15l we 
also plot the energies of bound states of an A atom with one, 
two or three B atoms. The vertical lines with arrows indicate 
the energies of few body bound states AB n . From left to 
right are shown the solid magenta line: energy of the tetramer 
AB3 Eq. (|47[) with three B atoms separated by distances riz = 
r23 = d and ri2 = V2d, dash-dotted black line: energy of a 
AB2 trimer Eo. (|46|) with two B atoms separated by a distance 
ri2 = d, dash-dotted blue line: energy of AB2 with two B 
atoms separated by a distance ri2 = V^d, black solid line: 
energy of a AB dimer Eq. (|45p . dashed blue line: energy of 
AB2 with two B atoms separated by a distance r\2 = V^d. 



In 3D, we found evidence for the occurrence of sev- 
eral energy mobility edges, for both positive and nega- 
tive matter wave energies E. For E > 0, we found a 
mobility edge for a positive effective scattering length of 
the order of the mean distance between scatterers. For 
a too large positive, or for a negative value of a c ff we 
found no evidence of any mobility edge, in agreement 
with the predictions in [14J. For a c g small and posi- 
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tive we found the occurrence of an energy gap between 
E = and an upper bound coinciding with a mean field 
effect E = g e sp where p if the density of scatterers, and 
5cff = 2irh 2 a c ff/m. Within this energy gap, there is a 
small localization length with no evidence of localized 
states [HI]. For E < 0, where the matter wave is bound 
inside the gas of scatterers (a case not explored in previ- 
ous works to our knowledge), we found localized bound 
states and evidence of two mobility edges separated by 
an energy interval where only extended bound states are 
present. We also found, for a e s small and negative, ex- 
tended bound states in the energy interval between the 
mean field shift E = pg c g and E = 0. To ascertain 
the presence of the mobility edges, and to identify the 
universality class of this disorder induced transition, a 
dedicated analysis should be undertaken using finite size 
scaling techniques to determine the critical exponent. For 
experimental purposes, large values of the effective scat- 
tering length may be obtained by using Feshbach reso- 
nances to control the free space scattering length aj9j. 
Recently a > lpm has been obtained with Li gases [331 ] . 

In 2D, no evidence of mobility edges is found for either 
positive or negative energies. Contrary to the 3D case, 
at negative energies, all bound states are localized. We 
identified regions in the (E/Eo, ln(a e ff/d)) plane where 
the localization length is small and the density of states 
is high. This happens for a c g of the order of the average 
distance between scatterers, for a wide interval of ener- 
gies. At high energies we find a rapid, exponential-like 
increase of the localization length with energy At nega- 
tive energies the localization length shows no such rapid 
increase of £. To our knowledge, large values of the 2D 
effective scattering length have not yet been observed. It 
is thus interesting to study the dependence of the 2D a c s 
with respect to the free space scattering length a, and 
the harmonic confinement in the microtrap. 
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FIG. 17: For the 2D system: For all the resonances with 
Er CB /E £ (0,2) and ln(a cff /d) £ (-3,3), and for all the 
bound states with E/Eo £ (—2,0) and ln(a e ff/d) £ (—3,3), 
the figures show an histogram giving the number of resonances 
(in (a) and (b)) and bound states (in (c) and (d)) per scat- 
terer (the number of scatterers is N — 7 x 10 3 ) as a function 
of: in (a) and (c), the square root of the participation sur- 
face defined in Eq. (|44|l : in (b) and (d), the r.m.s. size in real 
space defined by the 2D version of Eq. (|31[) . The bin size is 
~ 0.6d for (a), ~ 1.37d for (b), ~ 0.3d for (c), and ~ 1.37d for 
(d) . The parameters are the same as in Fig. 114b . The dashed 
vertical lines are the values Sp /d ~ 9.5 and a/d ~ 4.2 used 
in Fig. I14fc and Fig. I14H to select the bound states and reso- 
nances that are spatially localized and to filter out the bound 
states and resonances that are weakly localized. 



matrix dM(E)/dE is positive in 3D (respectively neg- 
ative in 2D), so that drrii/dE > in 3D (respectively 
drrn/dE < in 2D). 

Let us start with the 3D case. The idea is to show 
that, for all 1 < i, j < N, 



We acknowledge useful discussion with D. Dclande. 
Numerical calculations have been performed on the 
IFRAF cluster (the entire work needed several CPU 
months per core, on 96 Intel Xeon Quadcore proces- 
sors). M.A. acknowledges financial support from the 
ERC Project FERLODIM N.228177. D.H. acknowledges 
financial support from CNRS, UPMC and IFRAF during 
his stay in Paris. 



dE 



Mij(E) 



m (E — ho)'- 



(Al) 



where ho = — ^-A^ D . Since (E — h ) 2 is a positive oper- 
ator, the positivity of dM /dE readily follows. For i ^ j 
one has, from Eq.([8j), that 



Appendix A: Monotonic behavior of the eigenvalues 

of M(E) 

1,2 2 

We show that, for E = — -g^- < ^' tne eigenvalues 
rrii(E) of the real symmetric matrix M(E) are mono- 
tonic functions of the energy E. This results from the 
Hcllmann-Fcynman theorem, and from the fact that the 



2irh 2 



^E-ho 



(A2) 



since (E — /io)go( r ) = ^( r )- Taking the derivative with 
respect to E gives Eq. (|Al[) . It remains to check that 
Eq. (|Al[) also holds for i = j by a direct calculation. On 
one hand, Ma = — q so that dMa/dE = m/(h 2 q). 
On the other hand, introducing a closure relation in the 
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plane wave basis, one indeed finds 
2nh 2 , 
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for all 1 < i,j < N, with Iiq 
has from Eq.® that 



1^2° For i 7^ j, one 
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since (E 1 — h )gn(r) = 5(r). For i = j, one again per- 
forms a direct calculation. First, dMa/dE — d\wq/dE = 
1/(2E). Second, a closure relation in the plane wave basis 
indeed gives 



irh 2 



1 



'(E-ho) 2 

irh 2 r d 2 k 

rn 
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(2tt) 2 (E-H 2 k 2 /2m) 
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A simple consequence of the monotonic behavior of the 
eigenvalues to,, is that the total number of bound states, 
for a given realization of disorder, is given by the number 
of positive eigenvalues of M(E = 0) in 3D, and by the 
number of negative eigenvalues of M(E — ^ 0~) in 2D. 




E/E. 







FIG. 18: (Color online) For the 2D system: Width Y of the 
resonances as a function of the energy E and the logarithm 
of the effective scattering length ln(a e g/d) [see Eqs. (|24|) and 
(|27|l , applied to the eigenvalues mf (E) of the matrix Moo of 
Eq. (|4HI ], The physical parameters are the same as in Fig. 1141 
The plane (E, ln(a c ff/d)) is decomposed in rectangular bins 
of widths 8E and 8 ln(a c ff / d) . The color map (see bar on the 
right) is applied to the quantity log 10 h{T) /Eo where (r) is the 
mean value of T for the resonances within a given bin. The 
resonances are filtered in (a) over the participation volume as 
in Fig. 114b and in (b) over the r.m.s. size a as in Fig. I14H. In 
(a) one has 8E = 0.019£ and 8 ln(a cff /d) = 0.027, and in (b) 
one has SE = 0.011£ and 8\n{a eB /d) = 0.024. 



In 2D, the proof is quite similar. One simply has to 



Appendix B: Matter wave wavefunction in terms of 

A 

We restrict here to the three-dimensional case, the gen- 
eralisation to the two-dimensional case being straightfor- 
ward. Let us consider first a matter wave bound state 
with eigenenergy Eq < 0. Eq is a pole of the Green's 
function G(r, ro) defined in Eq.([2]), so that, according to 
Eq.©, an eigenvalue mo(E) of the matrix M(E) van- 
ishes for E = Eq. Since M(E) is real symmetric, its 
inverse has the spectral decomposition 



N-l 



1 



( E ) = E 7^Ej D *W D i( E ) ( B1 ) 



71=0 



where (£>j n) ) i<i<N is the orthonormal eigenvector of 
M(E) with real components and eigenvalue m n (E). In 

particular, J^iLii^^ ) 2 = 1. When this spectral decom- 
position is injected in Eq.©, together with the expansion 
mo(E) = (E - Ea)m'(Eo) + . . ., it leads to 



G(r,r ) 



E^Eo 



y>o(r)y>o(r ) 

E — En 



(B2) 



with 



(2^/m)V^ ,(o )( ^ o(r _ rj) . (B3) 



[m' (E ))V2 



As it appears from Eq. (|B2[) . ip is the wavefunction of the 
bound state of energy Eq. Note that m' (En) is positive 
according to Appendix [AJ 
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Let us consider now a a resonance z Tes of the system 
with a complex energy z rcs . Then z rcs is a pole of the 
analytic continuation of the Green's function from the 
upper half-plane to the lower half-plane, so that, accord- 
ing to the corresponding analytic continuation of Eq. (J9j) , 
an eigenvalue of the matrix M(z) vanishes for z = z rcs . 
Since M(z) is complex symmetric, if it is diagonalizable 
its inverse has the spectral decomposition 

N-l 

[M-%{z) = £ -£—D < f\z)Df\z) (B4) 

n=0 n ^ ' 

where (p\ )i<i<jv is the right eigenvector of M(z) with 
complex components and eigenvalue ra„(z), the corre- 
sponding left eigenvector is simply its complex conjugate, 
so that the normalisation condition is J^iLii-D^) 2 = 1. 



When this spectral decomposition is injected in the ana- 
lytic continuation of Eq.([9|), together with the expansion 
mo(z) = (z — z rcs )m' (z Tes ) + . . ., it leads to 

G(r,r ) ~ Mr)Mr ) (B5) 

Z—¥Z T ea Z Z res 

with 

= gS^ E^ 0) (^) g o(r-r,). (B6) 

In the limit of an infinitely extended disordered system 
-00 would correspond to the wavefunction of a localized 
state of energy z res . 
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